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Abstract 

A tight binding parameterization of the band structure, along with a mean 
field treatment of Hund, electron-electron, and electron-lattice couplings, is 
used to obtain the full optical conductivity tensor of LaMn03 as a function of 
temperature. We predict striking changes with temperature in the functional 
form and magnitude of the optical absorption. Comparison of our results to 
data will determine the Hund, electron-lattice, and electron-electron interac- 
tions. 
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I. INTRODUCTION 



In recent years the colossal magnetoresistance manganese perovskites have attracted 
much attention.! The low energy (hu> <4 eV) physics of these materials is believed to be 
governed by Mn e g electrons, which are coupled by a strong Hund's coupling, J H , to Mn 
t2 9 symmetry 'core spins' and also interact with each other and with lattice distortions. 
Three important and still unresolved issues are the value of Jh, the strength of the Jahn- 
Teller(JT) electron-phonon coupling A, and the magnitude of the effective e g -e g on-site 
Coulomb interaction U governing the low energy physics. The relative importance of the 
electron-electron and electron-phonon interactions is a particularly important open question. 

Hund's coupling Jh is believed to be large, because in good samples at low temperatures, 
observed magnetic moments are close to their saturation values.! However, Jh has not been 
directly measured. Estimates ranging from 2JnS c ~ 2 eV (Ref. to 3 eV (Ref. f4|) have been 
obtained from the analysis of the optical measurements of doped materials. However, these 
estimates depend upon proper identification of features in the observed optical conductivity, 
and this identification is at present disputed. 

The value of the effective Coulomb repulsion is also unclear. High energy photoemission 
spectroscopy! yields a bare repulsion [/bare ~ 8 eV, but also implies that the manganites 
are in the class of charge-transfer materials,!^ as are the high-T c cuprates. In charge-transfer 
materials the low energy [Twj <C L^are) electronic degrees of freedom are complicated linear 
combinations of transition metal and oxygen orbitals, and the effective repulsion for these 
orbitals may be much less than [/bare- F° r example, in high T c materials, [/bare ~ 10 eV 
but the U relevant to the low energy physics is rather less; perhaps (7^2 eV (Ref. 0). 
The proper value of the U relevant to the low energy {Two < 4 eV) physics of manganites is 
at present disputed. Some authors argue that electron-electron interactions are dominant; 
others that electron-phonon interactions are crucial. 

In this communication, we show how measurements of the magnitude, anisotropy and 
temperature dependence of the optical conductivity of the manganites' insulating 'parent 
compound' LaMnOs can help to resolve these issues. We obtain an effective Hamiltonian 
for the low energy degrees of freedom by fitting band theory calculations to a tight binding 
model and adding Coulomb interaction terms. We then determine its optical conductivity as 
a function of temperature, and show how the different features in the spectrum may be used 
to determine the interaction parameters. The advantage of the tight binding approach is that 
it compactly represents the low energy degrees of freedom, allows straitforward inclusion of 
interaction effects, and is computationally cheap enough that many parameter choices may 
be investigated. 

The material is a (0,0, tt) antiferromagnet at low temperatures, while a paramagnet at 
room temperature. The Hund's coupling therefore leads to a peak structure with a pro- 
nounced and strongly temperature-dependent anisotropy. Observation of this effect would 
provide a definitive determination of the Hund's coupling parameter. 

The Jahn- Teller and Coulomb interactions lead to similar effects in the absorption spec- 
trum; indeed the two interactions roughly add. The strength of the Jahn- Teller coupling 
may be determined from crystallography data or from band theory; interestingly the two 
determinations agree very well. The effect of the Jahn- Teller coupling on the optical conduc- 
tivity may essentially be calculated exactly. We compare the results of this calculation to 
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the optical conductivity data, and argue that differences are due to the Coulomb interaction. 
To estimate the strength of the Coulomb interaction we use a Hartree-Fock approximation 
and calculate the optical conductivity of the Jahn- Teller plus Coulomb system. The opti- 
cal conductivity of LaMnOs was previously calculated by Solovyev et a/.l and by Terakura 
et a/.i using the LDA method, but they did not consider the temperature dependence, 
or present information about the dependence of the observed conductivity on parameters. 
When we calculate using parameters tuned to give similar band gaps, we find the maximum 
conductivity is factor of 2-3 larger than the result in Ref. [8] and 4 times than in Ref. |9|. 
The difference in absorption strength is not understood at present, but is discussed in more 
detail in Section V, which shows that the error is unlikely to arise from the tight binding 
approach. 

The rest of this paper is organized as follows: Section II presents the Hamiltonian, band 
structure, and the parameter determination, Section III gives the formulation of the optical 
conductivity, calculation results, and the comparison with experiments, Section IV describes 
the effect of the Coulomb interaction, Section V discusses the accuracy of the approximation 
we have used, and Section VI is the conclusion. 



II. HAMILTONIAN AND BAND STRUCTURE 

At and below room temperature, LaMnOa exists in a distorted form of the ABO3 per- 
ovskite structure. According to band theory calculations0~@ the conduction band is derived 
mainly from Mn e g symmetry d-orbitals and is well separated from other bands. We find 
that the band structure appropriate to the ideal cubic ABO3 perovskite structure may be 
well represented by the following tight binding model. 

H KE + H li = — Y tfdl d 7 ,t h + H.c. -uY dl d 7nn . (1) 

i,(5, a, 6, a i.a.a 



Here i represents the coordinates of the Mn sites (which in the ideal structure are arranged in 
a simple cubic lattice), a, b represent the two degenerate Mn e g orbitals on a site, S(— ±x, y, z) 
labels the nearest neighbors of a Mn site, a denotes the spin state, and t^ b is the hopping 

amplitude between orbital a on site i and b on site i + S. We choose \ipi >= \3z 2 — r 2 > 
and \tp2 >= \ x<2 — > as the two linearly independent e g orbitals on a site. The hopping 
matrix tg b has a special form: for hopping along the z direction, it connects only the two 
|3z 2 — r 2 > states, thus 



t z = f_, = t n n . (2) 




The hopping matrices in the other bond directions are obtained by appropriate rotations 
and are 
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As noted in the introduction, there is a substantial high energy ([/"bare ~ 8 eV) photoemission 
evidence for strong on-site Coulomb interactions in the manganites, which places them in 
the class of 'charge transfer' materials. The relevance of the band theory calculation may 
therefore be questioned. We argue, however, that the effects of the interactions at the 
\ow(Hu < 4 eV) energies of interest may be determined by comparing the predictions of 
the band theory calculation to data; the results we present will allow this comparison to 
be made. The high T c superconductors provide an instructive example. These are also 
charge transfer insulators with a very large high energy on-site repulsion [/high (Ref. [7|,[13] ). 
The low energy excitations are complicated objects called Zhang-Rice singlets, but it has 
been established that the effective interaction relevant to the low energy theory is much less 
than Lehigh, and that band theory(albeit with a renormalized hopping) describes the electron 
dispersion well. We therefore suggest that band theory is an appropriate starting point in 
the manganite well. 

We now turn to the electron-lattice coupling. Below 800 K, LaMn03 exists in a distorted 
form of the ABO3 perovskite structure. The important distortion is a Jahn- Teller distortion 
in which Mn-0 bond lengths are changed in such a way that if we choose a Mn site as 
an origin, the distortion is even-parity and the sum of the lengths of the 6 Mn-0 bonds is 
unchanged. This distortion lifts the degeneracy of the e g levels on a site. To represent this 
we define u% as the displacement along the a direction of an oxygen ion located between Mn 

ions at % and i + a, and we define v " = tt? — u% „ . The Jahn- Teller distortion term may then 
be written as 

The experimentally observed distortion has two components: a Q2 type staggered distortion 
with wave vector (71*, 7r, 0), and a Q3 type uniform distortion. This distortion leads to a 
Jahn- Teller term of the form 





where w and v are the amplitudes of the staggered(Q2) and uniform(Q3) distortion respec- 
tively. 

We next consider the Hund's coupling. This leads to a term 

#Hund = J H^ • d la,*° a P d la,P' ( 7 ) 
i,a,a 

where S 1 * represents the t2 g core spin and a the Pauli matrix. At T = K, the magnetic 
structure is of a (0,0, n) antiferromagnet; leading to 

H„ U n d = JnS c E [(l " Mr) ia/l^ + + 4-Aa] • 

i,a 

The total Hamiltonian is the sum of the terms considered so far. 
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Htot — HkE + Hp + ifjT + -^Hund- 



(9) 



Fourier transforming into k space results in the expression for H tot given in the Appendix 
A. By diagonalizing this matrix, we can find the energy levels. 

At T = K, the unit cell is doubled twice, once by spin and once by orbital ordering. 
We have two orbital states for each of the two spin states, and four Mn sites per unit cell. 
Due to the symmetry between two spin states, we will have two-fold degeneracy for each 
level. Therefore, we have eight separate bands. The ground state is obtained by filling 
the energy levels below the chemical potential /i, which is determined to give the correct 
number of electrons per unit cell. We denote the energy levels by Ej(k) (j=l,2,...,8) in 
the order of the increasing energy. The band structure determined by Eq. (^[) is shown in 
Fig. 1 for the parameters which provide the best fit to the published band calculations.0 
Crudely speaking, the bands fall into 4 pairs, which may be understood by setting t Q = 
(as occurs at (tt/2, tt/2, tt/2)); in this case we have four separate energy levels on each 
site: which are i? lj2 = —Xy/v 2 + w 2 , E^^ = X\/v 2 + w 2 , E 5j6 = 2JuS c — Xy/v 2 + w 2 , and 
£7,8 = 2JhS c + X\/v 2 + w 2 . When we have non-zero t a , these levels split and become 
dispersive. In the low temperature (0,0,7r) antiferromagnetic structure, the Hund's coupling 
suppresses the z-directional hopping; the bands thus become more two-dimensional as JuS c 
increases. 

Band theory calculation is used to determine the three parameter values of our model 
Hamiltonian: £ > A, and JnS c . To find these, we fit our band structure calculation to the 
LDA band calculation for the JT distorted LaMn03 by Satpathy et alM at high symmetry 
points in reciprocal space, (7r,0,0), (0,0,0), (tt/2, tt/2, tt/2), and (7r,0,7r/2). The standard 
deviation is m 0.2 eV, and maximal error of 0.4 eV occurs at (tt/2, tt/2, tt/2) for the lower 
JT level of the upper Hund state, E 5t6 . The determined parameter values are t o =0.622 eV, 
A=1.38 eV/A, and 2JnS c =2A7 eV. The fitted band structure is shown in Fig. 1. These 
parameters fit the LDA band calculations for the JT distorted and buckled actual LaMn03 
structure published by Satpathy et a/.0 with a similar size of error. The above t Q and JhS c 
are similar to the values obtained by Myrasov et a/.0 from an LDA calculation for the ideal 
cubic structure. 

This parameter A may be independently determined by fitting the observed lattice 
distortioril3 to a simple model of localized electrons which are Jahn- Teller-coupled to a 
harmonic lattice as explained in Ref. [R| This reference shows that the amplitude of the 
observed distortion fixes the parameter X/(Kiao), where K\ is the spring constant between 
the nearest neighbor Mn-0 pairs and ao is the average Mn-0 distance. The following two 



equations are derived in Ref. 16 



°° ea = iz L t c ° s 2 ^ + + c ° s 2 (° 2 + ^ a )] ( io ) 

< = ^F [cos2 (0! + r) ~ cos 2 (9 2 + ^)] (11) 

where e a (a = x,y,z) is the uniform strain, it" is the staggered oxygen displacement, K 2 is 
the nearest neighbor Mn-Mn spring constant, ip x ' y > z = —tt/3, tt/3, 0, and Q\, 62 parameterize 
the orbital states on the two sublattices: \ip >= cos8\3z 2 — r 2 > +sin8\x 2 — y 2 >. (Eq. 
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([10|) corrects a factor of two error in Eq. (10) of Ref. |16| .) Using the relation 9 2 = tt — Q\ 
derived in Ref. 16], we get the following expression of X/(Kia ). 



Kia \ 



2K 2 /K X 



+ (\/3ao) 



(12) 



From Ellemans et a/.'s resultsjL3 e z 



0.0288, u x s =0. 141 A, and a =4.034 A. For < K 2 /K l < 
1.0, we obtain 0.0428 < X/(Kia ) < 0.0591. K\ is estimated from the frequency of the 
highest lying bond stretching mode measured in this material by Jung et a/..E3 The measured 
bond stretching mode has a peak at 70.3 meV. From the relation (hcu) 2 = 2Ki(m^ n + Wq 1 ), 
we obtain K\ =7.36 eV/A 2 . Using ao = 4.034 A, the range of A obtained is between 1.27 
eV/A and 1.76 eV/A, which includes the value obtained above. We can, in fact, determine 
the lower bound of K 2 /Ki from the structural transition temperature as explained in Ref. 
16|. In Ref. [16], the mean field estimation of the structural phase transition temperature was 



found to be T S MF = 3\ 2 K 2 /[2K 1 (K 1 + 2K 2 )\ = 750 K = 65 meV. Considering that mean 
field theory overestimates the transition temperature, we obtain 



2K X (K X + 2K 2 ) 



> 65 meV. 



(13) 



Combining the previous expression of \/{K\ao) versus K 2 /Ki, we can determine the range 
of K 2 jK\. The determined range is K 2 jK\ > 0.26, and gives A > 1.36 eV /A, which is 
remarkably close to the value obtained by band fitting. 



III. OPTICAL CONDUCTIVITY 



A. At T= K 



Optical conductivity per volume, a, can be found in the following way. The electromag- 
netic field couples to the electrons via the Peierls phase factor t$ — > t^exp (ieA- 5a /hJ. 
This implies the total current operator 



8A 



leao 
"2T 



E 



.,5,a,b,( 



r: Vi ieA ' 5a ° )S d- 
n 



+5ba 



H.c. 



(14) 



By expanding about A, we obtain 



J= J p + J d + 0(A 2 



(15) 



where 



- > ieao ^ t d 

2,(5, a, b, a 



H.c. 



a 2% ^— ' 5 % \ iaa i+oba J 

i,S,a,b,a 



(16) 
(17) 
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Since J(A) =< 0(A)\J p + J d \0(A) >= -iujUA + 0(A 2 ), where E = / dVa is the total 

conductivity, by expanding J (A) up to linear order in A, we can find a. Since Jd is already 
linear in A, it gives the following contribution in a. 

\v 1 fea \ 2 1 



V h ) 



where 



K 



AA 



N 



— < 01- Y tf (dt d u ? h „ + H.c.) 10 > . 



(18) 



(19) 



i,<5=±A,a,6.c 



Since J p does not contain A, we should use higher order perturbation to obtain linear term. 
Second order perturbation yields 



a 



1 



iuN Mn ao 



E 



< 0Ul\\n >< n|J o j0 > 



pi/ 1 



(20) 



where e is an infinitesimal introduced to make the expression well defined. The above 
procedure is a standard linear response theory, which is also used in Ref. 18, and explained 
in detail in Ref. [19. 

By transforming into energy eigenbasis, we can obtain explicit expression of K\\ and a p , 
which gives a. The results are 

2 



K xx = j—^ j R dk £ 2cosA; a fl a (£) ii 



Ej(k)<ii 



and 



cr 



pjaa 



iujo? (2ir) 3 Jr 



[ea /h)2 sin k a B a (k)jj> 

dk =; =; 

E 3 (k)<^EAk)>^ hu} ~ E A k ) + E A k ) + it 



where 



B a (k) = M(kYT a M(k), 



T 



T 



T 7 



(t x 








^ 





-t x 














t x 





V o 








-tj 


(ty 








\ 





— ty 



















V0 








-ty) 


/*, 








\ 





tz 











- 


-t g 





V o 








-tj 



(21) 



(22) 

(23) 
(24) 

(25) 

(26) 



7 



M(k) is the matrix diagonalizing H^^ in the Appendix A, a Q is the distance between Mn 
ions, and R is defined in the Appendix A. The factor of 2 comes from the two spin states. 
Therefore, the real part of a is 

1 ^ ' <l ° 2 sin k a B a (k)j 



Re[a aa ] = — sJtt^ / dk ^ 



Ej{k)<ii,E,,{k)>ii 



[huj - E r {k) + Ej(k)] 2 + e 2 

(27) 



The total oscillator strength in the conductivity of our Hamiltonian may be obtained as 
follows. 00 Optical conductivity a can be written as 



e 2 



where 



a Xv {u) = ( -5—K x J Xv - X xu ) /(-iw), (28) 



1 v <0|J p t A |n><n|J p ,|0> 
XXu N Mn al^ hu-{E n -E )+ie ' [ ' 

At large u, \\v oc \ju). Therefore, lim^^ a Xv {u) = j^K\\5\ u /(—iuj). Large ui limit of 
the Kramers- Kronig relation 



Im[a{u)} = -(2u/n)P f°° ds ^P^ (30) 

JO S z — LU Z 

yields the following sum rule: 



K^ = ^-l dujRe[a\\{uj)}. (31) 



e 2 7T Jo 

From crystallography studies in Ref. [15], w=0.488 A and u=0.174 A at T = K. Figures 
2(a), 2(c), and 2(e) show the T=0 K optical conductivities a xx and a zz calculated for three 
values of the coupling constant A with t Q and JnS c predicted by band theory and v and 
w from crystallographic data. Figure 2(a) shows a xx and a zz for the case 2\\/v 2 + w 2 < 
2JftS c . For a xx (solid line), we see a large peak at the Jahn- Teller splitting, corresponding 
to motion within one plane. Note the jump in absorption at the gap edge, characteristic 
two-dimensional feature, a weak feature at 2JuS c , corresponding to electron trajectories 
which overlap from one plane to the next, and an extremely weak feature at the sum of the 
Jahn- Teller and Hund's splitting. For a zz (dotted line), we see a very weak feature at the 
Jahn- Teller energy, corresponding to a small amplitude for an electron to tunnel through an 
intervening plane and land on a 'correctly oriented' core spin, a large peak at the Hund's 
energy, and another peak at the sum of Hund's and Jahn- Teller energies. The sharp peak at 
the Hund's energy in o zz originates from the essentially parallel bands seen in Fig. 1 between 
(tt/2, 7r/2, 7r/2) and (ir, 0, 7r/2). In LSDA band calculation in Ref. [11], these two bands are not 
exactly parallel, but deviates by 0.34 eV. Therefore, we expect this will cause the broadening 
of the peak by e=0.17 eV. The effect of this broadening is further discussed in Section V. 
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Figure 2(c) shows the case of comparable Hund's and Jahn- Teller couplings. We see that 
the structure becomes more complicated as the features overlap, and the band structure 
becomes less two-dimensional. The low energy shoulder starting from 1 eV originates from 
the transition between the opposite spin directions. Finally, Fig. 2(e) shows the case of the 
Jahn- Teller coupling greater than Hund's coupling. In this case the Hund's feature appears 
strongly for both a xx and a zz , whereas the Jahn- Teller feature is now almost completely 
absent in a zz . Figures 3(a), 3(c), and 3(e) show the results for a lower value of Jh-S'c, which 
show similar features as Figs. 2(a), 2(c), and 2(e). 

We have also studied the change of the optical conductivity caused by the change of the 
ordered orbital state. The orbital ordering angle, 9, defined in Section II, is related to the 
lattice distortion v and id by 

cos 29 = V (32) 

yv 2 + w 2 

For observed v and w at T = 0, 9 54°. We have varied the ratio between v and w without 
changing JT splitting 2X\/v 2 + w 2 . The results are shown in Fig. 4. When 9 = 7r/6, 
the orbital ordering is x 2 — z 2 /y 2 — z 2 type, and when 9 = n/3, 3x 2 — r 2 /3y 2 — r 2 type. 
When 9 = 7r/4, orbital state is between the two configurations, and v = 0, i.e., there is no 
uniform distortion. As orbital ordering changes from x 2 — z 2 /y 2 — z 2 to 3x 2 — r 2 /3y 2 — r 2 , total 
spectral weight of a zz has substantially decreased. If 9 is further varied toward 9=n/2, which 
corresponds to non-staggered x 2 — y 2 type orbital ordering, then the z-direction conduction 
becomes smaller. For a xx , the spectral weight has moved close to the lower edge without 
appreciable change in the total spectral weight. 



B. At T=300 K 

The model described above is appropriate to describe the optical conductivity when the 
spin and lattice degrees of freedom are perfectly ordered at T=0 K. If the temperature is 
non-zero, then the core spins will fluctuate from the ground state configuration. In the 
actual material, long range order is lost at Tn ~ 140 K. Therefore, by room temperature it 
is reasonable to assume that the core spins are completely disordered. In principle, similar 
considerations apply to the lattice degrees of freedom, but because room temperature is much 
lower than the structural transition temperature 800 K, we may neglect lattice fluctuation 
and assume a static JT distortion. 

To describe the system at Tn <C T <C 800 K, we develop the effective Hamiltonian in the 
following way. Instead of choosing spin basis along a fixed direction independent of sites, 

— * — * 

we choose ff on site i as the direction of e g electron parallel to the core spin on site i, and JJ- 
as its opposite direction. Therefore, the Hund coupling energy is 

Hn^ = 2JnS c Y.d\ a ^. (33) 

i,a 

and ifyr do not change their forms by the change of spin basis. We define the angle 
between the core spin directions on site % and on site i + 5 as 0*f + g, so that the kinetic energy 
is given by 
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H 



KE 



- L E tf 

sin 



COS f^) 4,^ + 008 ^^..* 



+ sin 



'6 



2 / "i<4 i+Sbty 

+ H.c. 



(34) 



This new representation of the Hamiltonian reduces to the previous H tot at T=0 K. At 
T 3> T N , will be completely random. Therefore, we argue it is appropriate to average 
the Hamiltonian by setting < cos(9^ + g/2) >= < sin(0^ + ^/2) >= 2/3, which gives the 
following effective kinetic energy term 



H 



eff 
KE 



~3^E l t d laf d i+m + d M d i+m 

i,6,a,b 

+ dt .d?,*,, + dt ..dr,?^ + H.c. 



(35) 



Transforming into k space yields the effective total Hamiltonian shown in the Appendix B. 

At room temperature, the core spin directions are fluctuating in time and space, which 
broadens the levels. To incorporate this physics, we introduce a phenomenological broaden- 
ing T, leading to the following expression of the optical conductivity. 



Re[a aa ] 
1 1 
ual (27r) 3 Js 

T/tv 



dk 



e\(eao/h)2 sin k a B' a {k) 



33' I 



r/vr 



/ / dEdE' y , 



[E-E 3 {k)Y + T* 
1 1 



dk 



E 



2r|(ea /ft)2sin k a B' a {k) 



33 



ual (2nr Js [ha, - E r {k) + E 3 (k)] 2 + W 

where 

B' a (k) = M'(k)XM'(k), 



(36) 
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ft x 


t x 





\ 


2 





-t x - 


-tx 




3 


tx 


o t x 









\0 


-t x - 


~t x 


/ 
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3 


ty 


ty 











~ty ^ 




/ 
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(37) 
(38) 

(39) 

(40) 
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M'{k) is the matrix diagonalizing if| in the Appendix B, S is defined in the Appendix B, 



and t Xj y )Z is the 2x2 matrix defined previously. T may be estimated from the root mean 
square fluctuation in the hopping; we find 



T « t a ^< cos 2 (#/2) > - < cos(fl/2) > 2 « (41) 

The general features of a xx and a zz at T ^> Tn are these: Because we have random spin 
directions along both x and z directions, both a xx and a zz show JT, Hund, and JT+Hund 
peaks. Due to the anisotropy of the lattice distortion, we still expect anisotropy in the peak 
intensity. The broadening due to random spin directions means the peaks become smoother 
than the T=0 K case. 

Optical conductivities calculated for room temperature are shown in Figs. 2(b), 2(d), 
and 2(f). For this calculation, we use the same A, t Q , and Jh<Sc as in Figs. 2(a), 2(c), and 
2(e), but we use the room temperature lattice parameters, which differ slightly from the 
K lattice parameters. We obtain w=0.417 A and v=0.155 A from Ref. [15]. As expected, 
the peaks are substantially broadened and indeed in Figs. 2(d) and 2(f) only two peaks 
are visible. Figures 3(b), 3(d), and 3(f) show similar results, obtained for the parameters 
used in Figs. 3(a), 3(c), and 3(e). The upturn of the optical conductivity at around zero 
frequency is an artifact of the crude consideration of the fluctuation in our model. 

We have calculated the variation of spectral weight with temperature and parameter 
values; the results obtained by Eqs. ( plD and (|3l"D are consistent, and shown in the Table I. 
It shows that at T=0 K, K xx sensitively decreases as A increases, but insensitive to JnS c , 
and K zz decreases as A or Jb_S c increases, which can be understood from the spin and 
lattice configuration at T=0 K. On the other hand, at T=300 K, both K xx and K zz have 
moderate dependence on both A and JnS c , which originates from the paramagnetic spin 
configuration. Temperature dependence of the total spectral weight is relatively weak, even 
though the spectral weight of each peak depends sensitively on temperature. Particularly, 
when the Jahn- Teller splitting is the much less than the Hund's splitting (as in Figs. 2(a) 
and 2(b)), it is possible to identify the lowest energy feature at both K and 300 K as 
arising from transitions between the parallel-spin but different Jahn- Teller states, and to 
determine the spectral weight in this feature. When we define 

K a ,jT = ^ / duja aa (uj), (42) 

7ie Z JJTpeak 

we find K^. K JT =0.271 eV, K%ffi =0.151 eV for Figs. 2(a) and 2(b), whose ratio is between 
1/2 and 2/3, as predicted in Ref. f|. The extra T=0 K JT spectral weight is pulled down 
from higher energy peaks as the spin disorder is decreased. It is also noteworthy that the 
peak shape is more asymmetric at T = K than at T = 300 K, due to the 2-d character. 

We now compare our results to data. This comparison is preliminary because the avail- 
able data disagree. Optical conductivity for polycrystalline LaMnOa was measured at room 
temperature by Jung et aL0 Because the crystal directions are random in polycrystalline 
samples, the observed quantity is <r av = 2a xx /3 + a zz /3, which we have also plotted in 
Fig. 5 (We assume here that the crystallite size is large.). Figure 5(a) shows the results 



for the parameter values determined from band fitting. The data in Ref. 21 exhibit two 
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main structures; a lower peak centered at 1.9 eV with maximum intensity 420 fi _1 cm _1 and 
K av ,jT = 2Kx,jt/3 + K z jt/3 ~ 0.115 eV, and a peak centered at around 4.5 eV with a much 
larger intensity. They attribute the 4.5 eV feature to e g - 2p transitions beyond the scope 
of our model, and assign the peak at 1.9 eV to JT-split e g -e g transitions within the parallel 
spin manifold, which we agree with. In this interpretation, the transitions to the reversed 
spin states are obscured by Mn-0 transition. Recently, room temperature optical reflectiv- 
ity spectra using a cleaved single crystal surface of Lai_ x Sr x Mn03 have been measured by 
Takenaka et a/.£3 Although it is referred to as a single crystal, we believe that the sample of 
LaMnOa is micro-twined. In their results, the Jahn- Teller peak appears at around 2.5 eV, 
having maximum intensity around 600 fi _1 cm~ 1 with width around 1.5 eV, corresponding 
to -Kav,jT « 0.141 eV, and the Mn-0 peak appears at 5 eV with maximum intensity 2800 
Q _1 cm _1 . Similar results were obtained by Okimoto et al.E Takenaka et al.'s results show 
a weak shoulder at 1.9 eV, which our model cannot explain. 

From Fig. 5(a), one sees that, if the Jahn- Teller interaction were the only important 
one, the observed lattice distortions would lead to a peak in cr av at 1.2 eV with maximum 
intensity « 1200 fl^cm -1 , width « 1.0 eV, and -Kav,jT = 0.223 eV. The maximum intensity 
or spectral weight is much larger than observed in either experiment and the peak position 
is lower. From Fig. 5(c), we see that the data of Ref. E2I may be approximately modeled 
by use of a stronger electron lattice coupling (or larger amplitude lattice distortion), which 
moves the peak to higher energy and reduces its spectral weight. The data in Ref. [H] are 
very difficult to reconcile with theory, because as can be seen from Fig. 3(b), one cannot 
simultaneously fit the peak amplitude and energy: choosing interaction parameters to fit 
the peak position leads to an amplitude, which is too large. The combination of peak energy 
and amplitude could only be explained if the actual hopping were rather smaller than the 
band theory value (say t Q ~ 0.4 eV rather than 0.6 eV). 

Further optical data would be very desirable (especially measurements at lower T). For 
the present we assume the Takenaka et al.'s data are correct, and consider their interpretation 
in more detail. We believe that the combination of band calculation and estimates from the 
crystallography data adequately fix the magnitude of the Jahn- Teller splitting. We therefore 
believe that the differences between the Takenaka et al. data and Fig. 5(a) are mainly due 
to the Coulomb interaction whose effects we study in the next section. 



IV. COULOMB INTERACTION 

We now add an on-site Hubbard-type Coulomb repulsion to our Hamiltonian. Because 
we have 2 orbital and 2 spin states on each site, we have in principle 6 different Coulomb 
repulsion terms, which may be generally written (n is the density operator) 

^Coulomb = EEE f/ («,/3),(a,6)% a , a %/3,6 ( 43 ) 

For simplicity we study this Hamiltonian in the Hartree-Fock approximation, which we 
believe is reasonably accurate for the simple quantities (peak position and spectral weight) 
important for our analysis. Corrections to the Hartree-Fock approximation are due to quan- 
tum fluctuations. We have compared Hartree-Fock to exact results for the case of strongest 
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fluctuations, namely the one- dimensional Hubbard model, and find agreement to within 
30% for spectral weig htS (peak positions are not available). We believe that in the case 
of present interest, the combination of three dimensionality, the large core spins, and the 
localization due to the electron-phonon interaction renders the Hartree-Fock approximation 
sufficiently accurate. 

In the approximation, one of the two density operators is replaced by its expectation 
value, which is determined self-consistently. The approximation explicitly breaks symmetry 
in spin and orbital space, so the issue of basis choice arises. We choose the orbital basis 
picked out by the observed lattice distortion and the spin basis picked out by the magnetic 
ordering. We refer the higher and lower lying orbital states by + and — respectively, and 
the spin states by ft and JJ. defined previously. The mean-field Hamiltonian may then be 
written as 

^Coulomb = E ^+4+%+ + C V4fr- < %- + U ^+ d k+ d ^+ + ^-4^-%-' ( 44 ) 
i 

where U^ t+ = {/^IR- < >+ < "4+ > + < >> etc -- 

Eq. ([P|) may be reorganized into a term proportional to the total e g density operator, 
which renormalizes the chemical potential and is of no interest, a term which couples to 
the total e g spin operator and changes the Hund's coupling, and terms which renormalize 
the local Jahn- Teller splitting in a manner which differs for electrons locally parallel and 
antiparallel to the core spin. Therefore, ifjT + i^Hund + ^couiomb can be cast into the following 
form. 

Hj T + #Hund + ^Coulomb = £ A^V^ 2 + W 2 {fl^ + - fl^_) + X'^V 2 + W 2 {fl^ + - hty_) 

i 

+2J^S C (% + + ), (45) 



where 



Uft+ — U-tt- 

W 



X 't = X + V7W=k > ( 46 ) 
2\ / v z + w 2 



A; = A+ y_ (47) 

2J' n S c = 2J U S C + l -(U^ + + Efy_ - U t+ - CV). (48) 

Because the hopping matrices are defined in \tpi > and \ip2 > orbital basis, we need to 
transform the basis. For this purpose, it turns out useful to define 

A; = X' av + 5X', (49) 
A; = X' av - 5X'. (50) 

If we transform into k space in > and \ip2 > basis, we obtain the Hamiltonian to calculate 
optical conductivity. Details are shown in the Appendix C along with the expression of the 
number operators. 

As the simplest case, we consider the case where C/( Q! , ( 8),(o,&) = U, independent of (a,/?) 
and (a, b). For this case, we use t Q , A, and JrS c obtained previously and determine the values 
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of U by calculating optical conductivity at T=300 K and comparing with experimental JT 
peak position. For the JT peak at 2.5 eV in Takenaka et aVs results, we obtain [7=1.6 eV. 
The obtained value of U is close to the difference of the experimental peak position and the 
calculated peak position for the U=0 case. The room temperature results are shown in Figs. 
6(b) and 6(d). As we increase the value of [7, the peak position shifts upwards by « [7, and 
the peak intensity decreases. For [7=1.6 eV, the calculated maximum peak intensity is 730 
fi -1 cm _1 , width 1.2 eV, and i^fv jt ~ 0.145 eV. Therefore, the calculated spectral weight 
is close to the observed spectral weight. With these determined values of [7, we calculate 
T=0 K results shown in Figs. 6(a) and 6(c), from which we can see the anisotropy and 
temperature dependence of optical conductivity. 

In the no-hopping case, the energy levels on each site are — X'^v 2 + w 2 , X'^y/v 2 + w 2 , 
2JftS c — X'^y/v 2 + w 2 , and 2J{ l S c + X'^y/v 2 + w 2 . Even though the finite hopping gives dis- 
persion to the energy levels, the peak positions are close to the energy differences between 
different levels. Therefore, the JT peak position is close to 2X'^\ / v 2 + w 2 , the Hund peak 
position is close to 2J' R S C + (XL — A^)\/u 2 + w 2 , and the JT+Hund peak position is close to 
2J{ 1 S C + (A^ + A^)%/u 2 + w 2 . To see the effect of different type of Coulomb interaction and 
other parameters, we use different values of X'», J^Sc and calculate optical conductivities 
without changing the value of XL. We observe that the JT peak position and spectral weight 
do not change very much. Therefore even if we use a more general Coulomb interaction, 
with [7( Qi , l 8),(a,6) dependent on the indices, as far as we fix the JT peak position by fixing XL, 
its spectral weight does not change very much. 

For the above model, we also calculate the variation of the level occupancies as a function 
of U. The results for the lowest lying orbital < ; _ > are shown in Fig. 7 for T=0 K and 
T=300 K, which shows that as U increases, the e g electrons are more likely to have spins 
parallel to the core spins and stay in the ground state of the local lattice distortion. The 
curves however show that the value of U required to fit Takenaka et aVs data does not 
change the ground state occupancy much. 



V. COMPARISON WITH LSDA CALCULATION OF OPTICAL 

CONDUCTIVITY 

Terakura et al. have calculated the T=0 K optical conductivity using optical matrix 
elements and energies obtained from their LSDA band calculation.^ Their conductivity is 
strikingly different from ours in two respects. First, the form is different: the sharp peaks 
we find are absent in their calculation. We suspect that the difference is due in large part 
to the 0.01 Ry ~ 0.14 eV level broadening employed in Ref. [£| Figure 8 shows the effect 
of introducing an artificial broadening, e, into our calculation; the result is to be compared 
with Fig. 2(a). This figure shows that as broadening is increased, the peaks diminish in 
amplitude and become more symmetrical (although there is always more asymmetry in our 
calculation than in Ref. ^). Particularly, the Hund peak in o zz is expected to have additional 
broadening of about 0.17 eV as mentioned in Section III. Therefore, if we include 0.01 Ry 
broadening assumed in Ref. |9|, the total broadening for the Hund peak in a zz will be e ~ 0.3 
eV, which explains the absence of this peak feature in Ref. |9|. Therefore, we believe the level 
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broadening along with the presence, in the calculation of Ref. of other bands, accounts 
for the difference in line shape. 

A far more serious discrepancy is the difference in spectral weights. The area under the 
lowest conductivity peak in Ref. || is about a factor of 4 smaller than in our calculation. This 
difference seems not to be caused by a trivial error: in our calculation, the kinetic energies 
obtained from Eq. fl2lD , from Eq. (^l|) (i.e., direct integration of a), or from the Hellman- 
Feynman theorem agree. According to Hellman-Feynman theorem, K = K xx + K yy + K zz 
can be found from the ground state energy E /N Mn by 

d(E /N Mn ) 

= ~ t0 dto ' (51) 

where to is the hopping parameter defined in Section II. At T=0 K, we have calculated E 
as the sum of energies of the filled bands, 



En 



N Mn (2tt; 



- [ dk\E 1 (k) + E 2 (k)] (52) 

J R 



and evaluated K using Eq. (pi]). The result obtained in this way is consistent with the 
results in Table I. 

We examine the size of the possible error due to the following two approximations we 
have made: First, we have assumed that the hopping between Mn ions, which originates 
from Mn-0 hopping, can be effectively represented without explicit consideration of O band. 
Second, we have used tight binding approximation. 

To study the effects of Mn-0 hybridization on the conductivity in the dominantly Mn 
bands, we consider a simple model of a 1-d Mn-0 chain along x direction. Each unit cell 
contains one Mn ion at position R^ n = n^ao with a d-orbital represented by d\, and one 
oxygen ion at position Rf = (rii + l/2)a with a p-orbital represented by p\. We consider 
a Mn-0 hopping of magnitude tMn-o an d choose sign to reflect the symmetry of the O 
p-orbital (the sign can be removed by change of k space origin). In addition, to model the 
Jahn- Teller distortion, we consider alternating periodic potential A on the Mn site. We 
represent the energy of the d-level relative to the p-level by V. For simplicity we assume 
spinless electrons. This can be represented by the following Hamiltonian. 

h = - tM "~° (4pi - pUi+i - 4+iPi + pUi + H - c ) 

+ J2^(- 1 Y d ldi-Vplp t . (53) 

We obtain the exact band structure and optical conductivity and compare this to the band 
structure and conductivity obtained from nearest-neighbor tight-binding fit to the two up- 
permost (Mn-dominant) bands. The difference turns out small. The effective tight binding 
Hamiltonian is 

^ = -yE { d Ui + i + 4^-1 + H-c) + E 4f(-l)^k (54) 
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By transforming into k space, we can find band structure for H and H c g. The band structure 
of H e s is simply given by 



E = ±^4$ B ca&k + ^f (55) 

For given tMn-o? A? and V, we can fit the band structure of H eS to that of H to determine 
£ eff and A eff . 

For tMn-o=2.0 eV, A=1.0 eV, we try V = 1.0 eV and 10 eV. The obtained d-bands are 
shown in Figs. 9(a) and 9(b) (7=1.0 eV and 10 eV, respectively) as solid lines, along with 
the best tight binding fit as dotted lines. Fitted parameter values are t e ff =0.88 eV, A e g 
=0.59 eV for 7=1.0 eV, and t cfr =0.35 eV, A cfr =0.93 eV for 7=10 eV. It shows that when 
7 =1.0 eV, the fitting has error of 10 % of d-band width (comparable to the error in the 
fits used in Section II - IV), but when 7=10 eV, the fitting has negligible error. For these 
two cases, we assume half filling, and calculate optical conductivities shown in Figs. 9(c) 
and 9(d). The insets show the integrated spectral weight, K(%uj) = (2/i 2 a /vre 2 ) f Q <j(u)du. 
In this calculation, we assumed the Mn-Mn distance is ao=4.034 A, and the cross sectional 
area perpendicular to the direction of the chain is a 2 ,. For 7 =10 eV, the two calculations 
give almost identical results. For 7=1.0 eV, the tight binding fit has about 25 % larger 
spectral weight. We therefore expect our Mn-only approximation yields errors ~ 25 %. 

Next, to estimate the error of the tight binding approximation (i.e., of the Peierls approx- 
imation to the optical matrix elements), we consider the following Kronig-Penney model. 



Id 2 °° 



y5(x - 2n - 1) + y5(x - 2n) 



(56) 



where Ai, A2 < 0, h=m e =e=l, and spinless electrons are assumed for simplicity. Once the 
values of Ai and A2 are given, we can find the band structure, the eigenstates ipk(x), and the 
optical conductivity, and compare this to the Peierls approximation. The band structure 
has two bands; at n=l/2, one is filled and the other is empty. Therefore, we can calculate 
the optical conductivity via 



1 1 
a = — - 



<2.k\£\l,k> 2 



u2ttJ-z [ u - E 2 (k) + E^k)]' + e 2 



(57) 



where 1 and 2 are the band indices. For Ai = —4, A2 = —5, we calculate the exact band 
structure, shown in Fig. 10(a) along with the best tight binding fit from Eq. fl5"4] ) (t e fj=1.33, 
A eff =0.6). It shows that the error is about 7 % of the total band width. Figure 10(b) 
shows the calculated optical conductivity for the exact K-P model and tight binding fit. 
The spectral weight of the tight binding fit is about 20 % larger than that of exact result, 
as shown in the inset. 

In our calculation of optical conductivity, our band fitting has error about 0.2 eV, which 
corresponds to about 5 % of total Mn e g f band width. So, we expect our approximation 
has similar size error in spectral weight as the two cases considered above. Therefore, we 
expect that our calculated optical conductivity may have overestimated spectral weight by 
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about 20 %, but not 400 %. Therefore, we believe that within this error, our approximations 
are valid. 

The relation between the kinetic energy and the optical spectral weight follows from 
the two assumptions of gauge invariance and reasonably localized d-electron wave functions. 
The success of the tight binding fit confirms this localized character. A tight binding param- 
eterization of the band structure has been used to study cr(u) in other correlated electron 
contexts,0@ and works well in high T c case. The apparent discrepancy found here for the 
manganite, is an important issue for future research. 

VI. CONCLUSION 

We have calculated the optical conductivity of LaMnC>3 and have shown that the available 
data are consistent with the band theory estimation for hopping parameter and e g level 
splitting. Our main prediction for the change in functional form and magnitude of a as T 
is decreased below Tn is contained in Fig. 6. Experimental determination of the Hund's 
coupling as well as final validation of the model must await definitive measurements of the 
magnitude and temperature dependence of a{ui). 

In conclusion, we comment on the implications of our results. First, we note that the 
estimate of the electron-lattice coupling derived from band theory is in good agreement with 
that derived directly from crystallography data. Second, we observe that the electron-lattice 
interaction by itself does not account for the magnitude of the gap or the spectral weight in 
the absorption spectrum. A Coulomb interaction U ~ 1.6 eV is also required. This value puts 
LaMnOs in the weak-intermediate coupling range: the Coulomb interaction is approximately 
40% of the full band width 6 t G ~ 3.6 eV. In the simple one-band Hubbard model, a Coulomb 
interaction of this size (relative to the band width) does not significantly affect properties 
(such as optical spectral weights) at reasonable dopings of order 0.2 or larger. The effects 
of this moderate Coulomb coupling on properties of models of doped manganites deserve 
further attention. Many authors have argued on the basis of photoemission datall that the 
Coulomb repulsion is large (5-10 eV); however as noted by the experimentalists themselves, 
because the manganites are charge-transfer rather than Mott-Hubbard materials (as are the 
high-Tc superconductors) the U measured in photoemission is not directly relevant to the 
low (Two < 4 eV) energy physics of interest here. 

Our data analysis focuses on robust features (peak positions and spectral weights) and 
is insensitive to the approximations we made. Uncertainties in the tight binding parame- 
terization of the band structure leads to an errors 0.2 eV in peak position, which is not 
important here; the consistency of the peak position and spectral weights leads us to believe 
the band theory estimates of t D are reasonably accurate. Uncertainties in the estimates 
of electron-phonon coupling A could change our estimated Coulomb repulsion by around 
0.2 eV. We note however that we have not included any excitonic effects arising from the 
first neighbor interactions such as those proposed by Maekawa et a/..@ At T=0 K, the 
two-dimensional character and general flatness of the bands suggest that these might be 
important and interesting to look for. 

The crucial prediction of the present model is the dramatic change in optical absorption 
with temperature. This change is a robust feature of the model, and comes from a dramatic 
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shift in spectral weight caused by ferromagnetic spin ordering, along with a very nearly two 
dimensional character of the bands at T — > 0, caused by the between-planes antiferromag- 
netism. Early data! reported only a weak temperature dependence of the optical absorption; 
if these data are reproduced, then our fundamental picture of the manganites based on e g 
electron with electron-lattice and electron-electron interactions must be modified. 

Finally, we note that a troubling discrepancy with band theory calculations exists. Fur- 
ther work is needed to find origin of the differences. 

We thank H. Drew, S. Louie, O. Myrasov, and M. Quijada for helpful discussions, and 
NSF-DMR-9705182 and the University of Maryland MRSEC for support. 



APPENDIX A: 



Without considering the Coulomb repulsion, the total Hamiltonian at T- 
by the following expression. 

H tot = Y $ r H rd t, 



where 
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w 



-Aw 



\ -\w J ' 

R = {k\\k x \ + \k y \ < rr and \k z \ < n/2}. 



(A12) 
(A13) 



APPENDIX B: 



Without considering the Coulomb repulsion, the effective total Hamiltonian at T=300 
K is given by the following expression. 



cff,300K 
tot 



fces 



where 



fc V l, fc.fr' 2,fc,fr' l,fc+(7r,7r,0),fr' 2,fc+(ir,7r,0),ft' 
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Mi, M 2 , y, jy and G are defined in Appendix A. 
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APPENDIX C: 



At T=0 K, the total mean field Hamiltonian with the Coulomb interaction, 

tot ak a> k Q,fe' 

a=T,l,fce-R 



(CI) 



has the same form as H tot in the Appendix A with Jh'S'c — ► Jh.Sc, A — > A^ and the following 
additional term. 



ttMF _ V"^ rf fiMF r 

^add,T=0K - 2^ a a ,k n a dd,a,k a a,V 



(C2) 



a=T,J.,fceii 
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<5y a sw a o o 

V^H/ a 5y a o o j 



(C3) 
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sv n 



SX' a v 
-SX'v 







-SX'w 



-SX'w 



5\'j = SX' , 

sx\ = -sx'. 



(C4) 

(C5) 

(C6) 
(C7) 



Similarly, for T=300 K the effective mean field total Hamiltonian with the Coulomb 
interaction, 



^MP,300K = E rfttf MF.SOOK^ 



(C8) 



fees' 



has the same form as if t e ot ' 300K in the Appendix B with JftS c — > Jh^cj ^ ~~ an d the 
following additional term. 



ttMF,300K 
-"add 



(C9) 



fees 
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if 
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With these Hamiltonians for given f/( a ,/3) i ( a ,b), we can repeatedly calculate < > until 
its value converges. At T=0 K, < h a>a > is given by the following expression. 



(2^)5 j R d% ^[ P ^Q^(k)\n + [AW^^L', (C13) 



where 
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and P a (fc) is the matrix diagonalizing H a £. 
AtT=300 K, 
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FIGURES 



FIG. 1. Fitted e g band structure of LaMn0 3 : to = 0.622 eV, 2 J H 5 C = 2.47 eV, A=1.38 eV/A, 
and fi=0A eV. (tt,0,0), (0,0,0), (7r/2,7r/2,vr/2), and (ir,0,ir/2) points correspond to M, T, R, and A 
points in Ref. 11 respectively. 

FIG. 2. Optical conductivities a xx (solid lines) and <r zz (dotted lines) for to = 0.622 eV, and 
2JuS c = 2.47 eV without Coulomb repulsion(U=0). 

FIG. 3. Optical conductivities a xx (solid lines) and a zz (dotted lines) for to = 0.622 eV, and 
2Jh<S c = 2.10 eV without Coulomb repulsion (U=0). 

FIG. 4. Optical conductivities o xx (solid lines) and o zz (dotted lines) for to = 0.622 eV, and 
2Jh<S'c = 2.47 eV, A = 1.38 eV/A, U=0 with different orbital ordering angle 9. 

FIG. 5. Average optical conductivities cr av for t = 0.622 eV and 2 J H 5 C = 2.47 eV at T=300 
K without Coulomb repulsion(U=0). 

FIG. 6. Optical conductivities for to = 0.622 eV, 2J H 5 C = 2.47 eV, A=1.38 eV/A, and U=1.59 

eV. 

FIG. 7. Occupancies of the lowest-lying orbital versus U for to = 0.622 eV, 2Jh5" c = 2.47 eV, 
and A=1.38 eV/A at T=0 K and 300 K. 

FIG. 8. Optical conductivities a xx (solid lines) and a zz (dotted lines) for to = 0.622 eV, and 
2J U S C = 2.47 eV, A = 1.38 eV/A, U=0 with different broadening e. 

FIG. 9. Mn d band structures [(a), (b)] and optical conductivities [(c), (d)] of Mn-0 chain in 
a model explicitly considering O p-level (solid lines) and in the best-fit effective Mn-Mn chain 
model (dotted lines) for i Mn _ o =2.0 eV, A=1.0 eV, and V = 1.0 eV, t eff =0.88 eV , A cfT = 0.59 
eV [(a),(c)], and V=W eV, t cS = 0.35 eV, A eff = 0.93 eV [(b), (d)]. The insets show the integrated 
spectral weight represented in terms of kinetic energy versus photon energy. 

FIG. 10. Band structures [(a)] and optical conductivities [(b)] for exact K-P model[Eq. (|56|)] 
with Ai = —4, A2 = —5 (solid line) and for the best tight binding fit (dotted line), t e g=1.33, 
A e g=0.6. The inset shows the integrated spectral weight represented in terms of kinetic energy 
versus photon energy. 
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TABLES 



TABLE I. Total spectral weight for io=0.622 eV without Coulomb interaction (U =0), expressed 



in terms of K x and K z . 



2J H S c (eV) 


A(eV/A) 


K° x K (eV) 


K° z K (eV) 


K 3 £ 00K (eV) 


Kf 0K (eV) 


2.10 


1.38 


0.295 


0.242 


0.290 


0.211 


2.10 


2.20 


0.236 


0.206 


0.248 


0.179 


2.10 


2.90 


0.196 


0.184 


0.219 


0.155 


2.47 


1.38 


0.294 


0.217 


0.277 


0.199 


2.47 


2.20 


0.235 


0.185 


0.236 


0.167 


2.47 


2.90 


0.195 


0.165 


0.207 


0.144 
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Figure 1, Ahn et . al 
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Figure 2, Ahn et. al 
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Figure 3, Ahn et. al 
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Fig. 4, K. Ahn et al. 
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Fig. 5, Ahn et. al 
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Figure 6, Ahn et. al 




Figure 7, Ahn et. al 
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Figure 8, K. Ahn et al. 
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